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ABSTRACT 

We propose a new diagnostic of the "dusty torus" model for Seyfert nuclei. Dust grains in the torus are 
heated by the nuclear continuum, and reradiate mostly in the mid-infrared wavelengths. From the torus 
geometry, it is predicted that the emission at A < 10 /xm has strong dependence on the viewing angle. 
Since the dependence is predicted to be insignificant at A > 10 /xm, we study the flux ratio between 3.5 
/im (L band) and 25 /im; i?(L,25) = log[(:^3.5^m 5*1/3 5^„)/(i^25/im 'S',y25^„J]. three different samples 
(optically selected, X-ray selected, and infrared selected samples) of Seyfert galaxies, the observed values 
of R{L^ 25) between type 1 Seyferts (Sis) and type 2 Seyferts (S2s) are found to be clearly separated; 
i?(i, 25) > —0.6 for Sis while R{L, 25) < —0.6 for S2s. This implies universality of their torus properties. 
With this result and the other observational characteristics, we investigate the most plausible torus model 
among those presented in Pier & Krolik (1992, 1993). 

Subject headings: galaxies: Seyfert - infrared: galaxies 



1. INTRODUCTION 

Dusty tori around active galactic nuclei (AGNs) play 

. an important role in the classification of Seyfert galax- 

' ies. (Antonucci & Miller 1985; see also Antonucci 1993 

, for a review). Seyfert galaxies observed from a face-on 

■ view of the torus are recognized as type 1 Seyferts (Sis) 
' while those observed from a edge-on view are recognized 
, as type 2 Seyferts (S2s). In this way, the dusty tori act 

■ as a material anisotropically obscuring the emission from 
\ their interior region. 

Dusty tori themselves are also important emitting 
' sources in AGNs. Dust grains within the torus absorb 
, high-energy photons from the central engine, and re-emit 
1 them in the mid-infrared (MIR) regime. Therefore, in- 
frared radiation from the dusty torus emission is useful in 
studying the physical properties of the tori in AGNs (e.g., 
Dopita et al. 1998 and references therein). Since the tori 
are quite optically thick, the MIR spectrum is predicted to 
have strong dependence on the viewing angle [Efstathiou 
& Rowan-Robinson 1990; Pier & Krohk 1992, 1993 (here- 
after PK92 and PK93, respectively); Granato & Danese 
1994; Granato, Danese, & Franceschini 1996, 1997]. When 
the torus is observed from a face-on view, its hot inner sur- 
face is seen and the emission at A < 10 /im is enhanced. 
When observed from a edge-on view, the emission at A < 
10 fim. is obscured and thus weakened. Heckman (1995) 
observed that the averaged ratio of Y-band (10 /im) flux 
to nonthermal radio flux is higher in Sis than in S2s (see 
also Giuricin, Mardirossian, & Mezzerre 1995). Heckman, 
Chambers, & Postman (1992) observed a similar enhance- 
ment in radio-loud quasars (i.e., type 1) with respect to 
radio galaxies (i.e., type 2). PK93 observed that flux ra- 
tios of L band (3.5 /im) to Y band in Sis are higher than 
those in S2s. Fadda et al. (1998) observed that the MIR 



spectrum is steeper (i.e., redder) in S2s than in Sis. How- 
ever, further details of the MIR emission from dusty tori 
are unknown. 

This paper proposes the flux ratio of L band to IRAS 
25 /im band as a new MIR diagnostic for the dusty torus 
model (§2). We compile the observational data from the 
literature (§3), compare the above ratios of Sis with those 
of S2s (§4), and discuss properties of the tori (§5). 

2. NEW MIR DIAGNOSTIC 

As stated above, the torus emission is expected to be 
more anisotropic at A < 10 fim. than at A > 20 /xm be- 
cause the visibility of the inner wall of the torus is highly 
viewing angle dependent. Therefore, it is of interest to 
compare Sis with S2s in a flux ratio between A < 10 /im 
and A > 20 jixn. Since the IRAS photometric data are 
available for most of the nearby Seyfert galaxies (Moshir 
et al. 1992), we adopt the flux ratio between L band and 
IRAS 25 /im band, 

R{L,2b) = log[(l/3.5pm5'i,3 5^„)/(l'25Mm5'i.2,^„)]. 

The basic concept of our new MIR diagnostic is schemat- 
ically shown in Figure 1. Since the viewing angle depen- 
dence is more significant at 3.5 /zm, S2s are expected to 
have lower values of i?(L, 25) than Sis. 

Here we note that PK 93 used the flux ratio between L 
and Y bands, R{L,N), to compare Sis with S2s. How- 
ever, the Y-band flux is affected by a silicate line at 9.7 
/im. When the torus is observed from a edge-on view, the 
silicate line is seen as an absorption line, and thus both 
the fluxes in L and Y bands are weakened. The difference 
between Sis and S2s in R{L, Y) is thereby expected not 
to be as prominent as that in R{L, 25). 



1 



NEW MID-INFRARED DIAGNOSTIC OF THE DUSTY TORUS MODEL FOR AGN 



> 



-I — I 1 1 1 III 




0.1 1 3.5 10 25 100 

Wavelength (fim) 




Fig. 1. — Basic concept of our MIR diagnostic. The left panel shows typical spectra of the torus emission for Sis (upper) and for S2s 
(lower). The right panel shows how the 3.5 jxm to 25 ixm. flux ratio yields the viewing angle toward the torus. 



3. DATA SAMPLE 

To perform a statistical analysis with the MIR diag- 
nostic defined in the previous section, we have compiled 
photometric data in L, N , and IRAS 25 fxm bands from 
the literature (e.g., Ward et al. 1987; Roche et al. 1991; 
Moshir et al. 1992; PK93). Since radiation form AGNs 
is anisotropic in most of the energy bands, it is difficult 
to construct a statistically complete sample. We instead 
adopt three samples chosen by different selection crite- 
ria. The first sample consists of the CfA Seyfert galax- 
ies^ (Huchra & Burg 1992), which provide a well defined 
collection of objects limited by the B magnitude of their 
host galaxies. The second sample is the one limited by 
the hard X-ray flux from 2 to 10 keV (Ward et al. 1987). 
Since hard X-rays arise from the central engine itself and 
are not affected seriously by dust grains, this sample is ex- 
pected to be fair at least for Sis. The third sample is taken 
from Roche ct al. (1991). This sample is not complete but 
composed of A^-band bright objects. For each object in 
this sample, Roche et al. (1991) observed an emission fea- 
ture at 11.3 /im, which allows us to examine the presence 
or absence of any circumnuclear star formation activities. 
Our CfA, Ward, and Roche samples contain 18 Sis and 6 
S2s, 20 Sis and 4 S2s, and 11 Sis and 11 S2s, respectively 
Some objects are included in more than one sample. In 
total, there are 31 Sis and 14 S2s. Their basic data are 
summarized in Table 1. 

Besides the dusty torus, several sources in Seyfert galax- 
ies contribute to the observed MIR fluxes (see below). 
We exclude galaxies where the contamination with such 
sources appears to be significant. These galaxies arc indi- 
cated in Table 1. The resultant final samples consist of 27 
Sis and 5 S2s. 

The IRAS 25 /im measurements were made with an 
aperture which is large enough to cover the entire galaxy 
(0f75 X 4.'6; Neugebauer et al. 1984). There could be con- 
tamination with the disk of the host galaxy. To find galax- 
ies where the disk emission dominates over the torus emis- 
sion, we use the compactness parameter CP at 10 /im (Dc- 
vereux 1987), CP = fee x S^^/S^^^^^. Here 5,^^, is the N 

^The preliminary analysis based on the CfA Seyfert galaxies was 



band flux, 8^12^^ is the IRAS 12 /im flux, and fee is the 
color correction factor, fee ~ 0.12 S',yj2^„^/5,y25^,^ + 1.04. 
This compactness parameter gives an estimate on the ra- 
tio of the small-beam flux to the entire flux at 10 /im. The 
CP value of each galaxy is given in Table 1. Some objects 
exhibit CP > 1. This is due to uncertainties in the mea- 
surement or time variation of the nuclear flux. In such 
cases, we give CP = 1. The MIR fluxes of galaxies with 
CP < 0.5 are likely to be dominated by the disk emission. 
These galaxies are not used in our following analysis. 

On the other hand, the L- and A'^-band data given in 
Table 1 were obtained with small apertures {(f) ^ 5" 10"). 
In these data, the contamination with the host galaxy is 
unlikely to be important. From iiT-band images of Seyfert 
galaxies, Kotilainen ct al. (1992) estimated the average 
light contribution from the host galaxy as 32 %. Zitelli et 
al. (1993) found that L-band images of Seyfert galaxies are 
more centrally concentrated than the AT-band ones. Hence 
the contribution from the host galaxy to the L-band flux 
is less than ^ 30 %. 

Seyfert galaxies often exhibit circumnuclear starburst 
activities, which could affect the MIR emission (see Keto 
et al. 1992 for the case of NGC 7469). Such objects are ex- 
cluded from our analysis. As a signature of the starburst 
activity, we use emission features in the 8 13 /im regime 
(Roche et al. 1991). They are due to transient heating of 
polycyclic aromatic hydrocarbon molecules (PAHs) by UV 
photons from OB stars. Since PAHs are destroyed by X- 
rays, PAH features are absent in genuine AGNs (see Voit 
1992 and references therein). 

We also exclude narrow- line X-ray galaxies (NLXGs), 
i.e., S2s with strong hard X-ray emission (Shuder 1980; 
Veron et al. 1980; Ward et al. 1987). The central engine 
of these galaxies is believed to be hidden not by a dusty 
tori but by the disk of the host galaxy. Most of NLXGs 
are actually edge-on galaxies (see Ulvestad & Wilson 1984 
and Keel 1980 for the cases of NGC 2992 and NGC 5506). 
Furthermore, Glass et al. (1981) reported that R{L, N) 
values of NLXGs are similar to those of Sis rather than 
those of S2s. 

reported in Murayama, Mouri, & Taniguchi (1997). 
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Table 1 

Infrared Properties of the Sample Seyfert Galaxies 



Name'' 


L 


N 


^•^12 urn 

(Jy) 


(Jy) 


'S'l'60 Mm 

(Jy) 


'S'l^lOO om 

(JyJ 


CP^ 


R{L, 2hf 


Sample'' 


Rejecf* 




(Jy) 


(Jy) 










Type 1 Seyfert 


2237+07 


0.016 


0.082 


0.1400 


0.3920 


0.900 


1.270 


0.81 


-0.54 


C 




3A 0557-383 


0.113 


0.347 


0.5289 


0.6846 


0.3223 


< 0.5589 


0.78 


0.07 


W R 




3C 120 


0.090 


0.220 


0.2860 


0.6350 


1.283 


2.786 


1.00 


0.00 


W R 




Akn 120 


0.0740 


0.139 


0.3191 


0.4099 


0.643 


1.084 


0.52 


0.11 


W 




ESO 141-G55 


0.0640 


0.184 


0.2420 


0.3522 


0.5741 


< 1.481 


0.92 


0.11 


W 




I Zw 1 


0.110 


0.390 


0.5118 


1.211 


2.243 


2.643 


1 


-0.19 


C R 




IC 4329A 


0.227 


0.760 


1.082 


2.213 


2.030 


1.661 


0.90 


-0.14 


W R 




MCG -6-30-15 


0.0877 


0.286 


0.3803 


0.8088 


1.087 


1.096 


0.97 


-0.12 


W 




MCG 8-11-11 


0.075 


0.296 


0.6394 


1.948 


3.005 


4.235 


0.65 


-0.57 


W R 




Mrk 79 


0.0543 


0.200 


0.3062 


0.7625 


1.503 


2.363 


0.87 


-0.30 


W 




Mrk 231 


0.360 


1.420 


1.872 


8.662 


31.99 


30.29 


1 


-0.53 


C R 




Mrk 279* 


0.032 


0.076 


0.1990 


0.2890 


1.200 


1.970 


0.46 


-0.10 


C 


CP < 0.5 


Mrk 335 


0.123 


0.210 


0.3021 


0.3777 


0.3433 


< 0.5673 


0.83 


0.37 


c 




Mrk 509 


0.113 


0.220 


0.3158 


0.7018 


1.364 


1.521 


0.91 


0.06 


W R 




Mrk 530 


0.025 


0.077 


0.1800 


0.1910 


0.8560 


2.140 


0.50 


-0.03 


c 




Mrk 590 


0.0469 


0.169 


0.1917 


0.2214 


0.4893 


1.457 


1 


0.18 


c w 




Mrk 766 


0.0559 


0.288 


0.3855 


1.295 


4.026 


4.658 


1 


-0.51 


c 




Mrk 817 


0.0619 


0.357 


0.3350 


1.175 


2.118 


2.268 


1 


-0.42 


c 




Mrk 841 


0.0407 


0.1445 


0.1924 


0.4726 


0.4593 


< 0.6176 


1 


-0.21 


c w 




Mrk 1040 


0.105 


0.294 


0.6104 


1.315 


2.555 


4.551 


0.63 


-0.25 


w 




NGC 3227 


0.0783 


0.263 


0.6671 


1.764 


7.825 


17.59 


0.54 


-0.50 


c w 




NGC 3516 


0.117 


0.239 


0.4258 


0.8937 


1.758 


2.259 


0.73 


-0.03 


c 




NGC 3783 


0.131 


0.400 


0.8396 


2.492 


3.257 


4.899 


0.67 


-0.43 


W R 




NGC 4051* 


0.077 


0.297 


0.8554 


1.590 


7.131 


23.92 


0.44 


-0.46 


C W R 


CP < 0.5 


NGC 4151 


0.344 


1.400 


2.080 


4.600 


6.720 


8.600 


0.88 


-0.27 


C W R 




NGC 4593 


0.081 


0.182 


0.3441 


0.8089 


3.052 


5.947 


0.70 


-0.15 


W 




NGC 5033* 


0.0487 


0.031 


0.9452 


1.148 


13.80 


43.85 


0.04 


-0.52 


c 


CP < 0.5 


NGC 5273 


0.016 


0.134 


0.1340 


0.2420 


0.9910 


2.030 


1 


-0.33 


c 




NGC 5548 


0.0986 


0.210 


0.4006 


0.7690 


1.073 


1.614 


0.67 


-0.04 


c w 




NGC 7213 


0.115 


0.261 


0.6063 


0.7421 


2.666 


8.177 


0.51 


0.04 


w 




NGC 7469* 


0.1594 


0.600 


1.348 


5.789 


25.87 


34.90 


0.69 


-0.71 


C W R 


PAH 


Type 2 Seyfert 


Circinus* 


0.701 


6.00 


19.58 


71.29 


248.7 


315.9 


0.45 


-1.16 


R 


CP < 0.5 


Mrk 266 


0.007 


0.306 


0.2307 


0.9765 


7.432 


11.07 


1 


-1.29 


C 




Mrk 348 


0.039 


0.300 


0.3080 


0.8347 


1.290 


1.549 


1 


-0.48 


R 




Mrk 533* 


0.046 


0.217 


0.6724 


1.896 


5.588 


8.146 


0.44 


-0.76 


C R 


CP < 0.5 


NGC 1068 


1.920 


18.0 


39.7 


85.04 


176.2 


224.0 


0.59 


-0.79 


C R 




NGC 1275 


0.078 


0.674 


1.069 


3.539 


7.146 


6.981 


0.91 


-0.81 


R 




NGC 2110* 


0.047 


0.198 


0.3488 


0.8397 


4.129 


5.676 


0.75 


-0.40 


R 


NLXG 


NGC 2992* 


0.057 


0.249 


0.594 


1.422 


6.941 


14.44 


0.56 


-0.55 


W R 


PAH, NLXG 


NGC 3079* 


0.073 


0.091 


1.523 


2.272 


44.5 


89.2 


0.07 


-0.64 


C 


CP < 0.5 


NGC 4388 


0.074 


0.404 


0.9964 


3.463 


10.24 


18.10 


0.59 


-0.82 


C R 




NGC 5506* 


0.313 


0.643 


1.282 


3.638 


8.409 


8.886 


0.69 


-0.22 


W R 


NLXG 


NGC 5929* 


0.0073 


0.0186 


0.360 


1.570 


9.450 


12.00 


0.08 


-1.48 


C 


CP < 0.5 


NGC 7172* 


0.132 


0.141 


0.4374 


0.7612 


5.712 


12.29 


0.40 


0.09 


W R 


CP < 0.5, NLXG 


NGC 7582* 


0.201 


0.877 


1.620 


6.436 


49.10 


72.92 


0.82 


-0.66 


W R 


PAH, NLXG 



'^Galaxies shown with asterisks are excluded from the analysis because of possible contamination to the infrared emission 
''Compactness parameter defined in text 

^i?(i,25)=l0g[(l/3.5 )/(!^25,^mS'^25Mm)] 

'^C: the CfA sample; W: the Ward sample; R: the Roche sample 
Reasons for excluding from the analysis 
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4. RESULTS 

Figure 2 shows frequency distributions of R{L, 25) of 
Sis and S2s separately for the CfA sample (a), the Ward 
sample (6), the Roche sample (c), and the total sample (d). 
The galaxies excluded in the previous section are shown 
by white bars. Only the galaxies shown by black bars are 
used in the following analysis. 

All of the Sis have R{L, 25) > -0.6 while most of the 
S2s have R{L, 25) < —0.6. The S2 which lies exceptionally 
at R{L, 25) > -0.6 is Mrk 348. This galaxy exhibits no 
silicate absorption feature at 9.7 /um (Roche et al. 1991). 
Since the silicate absorption is a common property of S2s, 
Mrk 348 is considered to be in a face-on view like usual 
Sis. The absence of the broad-line region in this galaxy 
could result from, e.g., obscuration of the central region 
by a small cloud. There is no significant difference in the 
distributions of Sis and S2s among the three samples in 
Figures 2a c. thus our samples is probably free of large ori- 
entation bias. If we apply the Kolmogrov-Smirnov (KS) 
test, the probability that the observed distributions of Sis 
and S2s originate in the same underlying population turns 
out to be 0.275 %. When galaxies shown by white bars 
are included, the distribution of S2s is different among the 
three samples. This difference is likely to come from the 
different sampling criteria. 

Figure 3 compares the R{L, 25) ratio with the nu- 
clear absolute B magnitude [Mb (nucleus)] for Sis. The 
Mb (nucleus) values are taken from Kotilainen, Ward, & 
Williger (1993) and Granato et al. (1993), and are used as 
a measure of the luminosity of the central engine. If the 
observed R{L, 25) ratio depends on the intrinsic huninos- 
ity of the central engine rather than the Seyfert type, there 
would be a certain relationship between R{L, 25) and Mb 
(nucleus). Since no clear correlation is seen in Figure 3, 
we conclude that the difference in the intrinsic nuclear lu- 
minosity does not affect the observed value of R{L, 25) in 
Sis. This conclusion is applicable to S2s because Sis and 
S2s are likely to have the same torus properties. 

Finally, wc show the frequency distributions of R{L, N) 
in Figure 4. The KS probability that the underlying popu- 
lations of Sis and S2s are the same is 0.0404 %. Although 
this value is smaller than that for the R{L, 25) ratio, the 
separation between Sis and S2s in R{L, N) (Figure 4) is 
less clear than that in R{L, 25) (Figure 2). This is because 
the A'^-band emission is affected by the silicate feature at 
9.7 yum. 

5. DISCUSSION 

The R{L, 25) values for the Sis are clearly separated 
from those for the S2s at the critical value of —0.6. This 
limits the extent to which the dusty torus can vary among 
Seyfert galaxies; such variations would add "noise" and 
cause overlap between the two types. Hereafter, we com- 
pare our results with theoretical torus models of PK92 and 
PK93, and investigate the model which agrees best with 
the observations. 

Figure 5 shows the geometrical configuration assumed 
in PK92 and PK93. The torus surrounds cylindrically 
around the central engine and the broad-line region. The 
semi-opening angle ^open is given by the inner radius a 
and the height h of the torus; ^open = tan~^{2a/h). The 
viewing angle i is defined as an angle between the rotation 



axis of the torus and the line of sight. The critical view- 
ing angle Zcr is defined such that the broad-line region is 
visible at i < i„- Since the actual torus should be clumpy 
and not have a sharp edge, we expect icr > ^open- 

The torus emission is parameterized by the three quan- 
tities in the models of PK92 and PK93; 1) T: the effective 
temperature of the inner wall of the torus, 2) a/h: the 
inner aspect ratio, and 3) Tr and r^: the radial and verti- 
cal Thomson optical depths. In the upper panel of Figure 
6, we show the theoretical i?(L, 25) values as a function of 
the viewing angle i for six dusty torus models of PK92 and 
PK93. In the lower panel, the observed R{L, 25) values are 
shown separately for Sis and S2s. For each of the models, 
the observed critical ratio, i?(L, 25) = —0.6, yields a crit- 
ical viewing angle. The results together with the model 
parameters are given in Table 2. 

The derived critical viewing angle ranges from 46° to 87° 
(Models 1, 4, 5, and 6). Since Models 2 and 3 do not give 
the critical viewing angle, these models are not appropri- 
ate for dusty tori in Seyfert galaxies. To proceed further, 
we have to compare the results with other observational 
properties of Seyfert galaxies (see below). 

Narrow-line regions of S2s often exhibit conical mor- 
phologies, which are due to shadowing of the nuclear ion- 
izing continuum by the torus. The observed semi-opening 
angle of the cone ^open(NLR) is thereby equal to the semi- 
opening angle of the torus ^open- Table 3 summarizes sta- 
tistical results from observations of conical narrow-line re- 
gions (Pogge 1989: Wilson & Tsvetanov 1994; Schmitt & 
Kinney 1996). These results indicate 6'opcn(NLR) ^ 30°. 
On the other hand, Model 1 has ^open = 11°. Thus this 
model is not appropriate for dusty tori in Seyfert galaxies. 

The critical viewing angle i can be estimated from the 
number statistics of Sis and S2s if we observe Seyfert nu- 
clei from random orientations on the statistical ground, 

^(Sl) . . . ^ 

where A(S1) and A^(S2) are the observed numbers of Sis 
and S2s, respectively (Miller & Goodrich 1990). Table 
4 summarizes the results for three different surveys of 
Seyfert galaxies (Osterbrock & Shaw 1988; Salzer 1989; 
Huchra & Burg 1992). The derived critical viewing angles 
ranges from 27° to 46°. Since Models 4 and 5 give too 
large critical viewing angles («cr > 80°), they are not ap- 
propriate for dusty tori in Seyfert galaxies. Consequently, 
among the six models of PK92 and PK93, Model 6 with 
^open = 31° and icr = 46° is the best torus model. 

The R{L, 25) values of the S2s lie between -1.48 and 
—0.6, which correspond to the viewing angles between 86° 
and 46°. On the other hand, i?(L,25) values of the Sis 
lie between —0.6 and 0.37. This range is not explained 
by Model 6. The locus of Model 6 in Figure 6 is drawn 
down only to i = 41°. Since the R{L,25) value at the 
smaller viewing angle is expected to be nearly constant, 
it would be impossible to reproduce the R{L, 25) ratio as 
high as 0.37. One possibility that explains this higher ra- 
tio may be the 3 /im bump often seen in type 1 AGNs. 
This bump may be attributed to thermal emission from 
hot dust grains with T ~ 1300 K (e.g., PK93). Because 
the models of PK92 and PK93 assumed that the host dust 
component is an additional source to the torus, R{L, 25) 
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log (V3.5^m'Sv,,n,yV25nm'Sv,5„„,) log (Vs.Sjmi^v, ,n„/V25^m5'v,5,„„) 

Fig. 2. — Histogram of R{L,25) for the CfA Seyferts (a), the sample of Ward et al. (1987) (6), the sample of Roche et al. (1991) (c), and 
the total sample (d). Galaxies shown by white bars are likely to suflfer from contamination and are not used in our analysis (see text and 
Table 1). 
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Fig. 3. — Diagram of R{L, 25) vs. absolute B magnitude of nuclei for Sis. Filled circles represent galaxies used in our analysis. Two 
open circles are Mrk 279 and NGC 7469, which are included in the sample but not used in our analysis (see text and Table 1). Nuclear B 
magnitudes were taken from Kotilaiiicn, Ward, & Willigcr (1993) and Granato et al. (1993) which measured magnitudes of unresolved nuclear 
components deconvolved from two-dimensional CCD images of Seyfert galaxies. The Hubble constant of 75 km Mpc~^ is assumed to 
obtain the absolute magnitudes. 



Table 2 

The Dusty Torus Models and the Derived Critical Viewing Angles 



Model 


rpa. 








e 

^torus 


i f 

f'CT 




(K) 






n 


n 


1 


1000 


0.1 


1 


1 


11 


44 


2 


1000 


0.3 


0.1 


0.1 


31 




3 


1000 


0.3 


1 


0.1 


31 




4 


1000 


0.3 


1 


1 


31 


87 


5 


800 


0.3 


1 


1 


31 


82 


6 


500 


0.3 


1 


1 


31 


46 



'^Effective temperature of the torus 

''Inner aspect ratio of the torus 

'^Radial Thomson optical depth of the torus 

'^Vertical Thomson optical depth of the torus 

''Semi-opening angle derived from the inner aspect ratio 

^Expected viewing angle of the torus 
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Fig. 4. — Frequency distributions of R{L,N) for the total sample. 



Observer 




Fig. 5. — Schematic illustration of the geometrical configuration of the dusty torus model. 
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Table 3 

Semi-Opening Angles of Ionization Cones of Seyfert Nuclei 



Reference 


A^Seyfert^ 


6»opcn(NLR) 


Poggc 1989 


4 


26° ± 11° 


Wilson and Tsvctanov 1994 


11 


32° ± 8° 


Schmitt and Kinney 1996 


12 


29° ± 9° 



'^Number of the observed Seyfert galaxies 



Table 4 

Semi-Opening Angles Derived from the Statistics of Type 1 and Type2 Seyfert Galaxies 



Reference 






fi/2'' 


icr(stat) 


Osterbrock and Shaw 1988 


6 


9 


0.125 


27° 


Salzer 1989 


9 


7 


0.20 


34° 


Huchra and Burg 1992 


25 


23 


0.435 


46° 



'^Number of type 1 Seyfert galaxies 
'^Number of type 2 Seyfert galaxies 

"^Number ratio between type 1 and type 2 Seyfert galaxies corrected for the completeness of the survey 



is underpredicted at small inclination angles. Although 
it is controversial whether this host dust component is a 
separate component from the torus or the inner surface of 

the torus, those hot dust grains lie close to the central en- 
gine in either case. Since their emission is important only 
when the central region is clearly visible, the 3 fim bump 
is negligible in S2s. Therefore the critical R{L, 25) and icr 
are not affected by the treatment of the 3 ^m bump. 

We have examined only the small and sparse sets of the 
model parameters presented by PK92 and PK93. Further 
analyses with the larger and denser parameter sets are re- 
quired to understand the torus properties in more detail. 
Nevertheless, our most important result, which has been 
obtained firstly with the new MIR diagnostic, is the clear 
separation in the R{L, 25) ratio between Sis and S2s. This 
strongly suggests that the torus properties do not vary 
among Seyfert galaxies. The effective temperature of the 
inner wall may be universal as a result of that the inner 
wall is formed by balancing the rate of dust destruction 
with the rate at which the torus clouds drift inward (Kro- 
lik & Begelman 1988; PK92). We suspect that there are 
also certain mechanisms confining the vertical structure of 
the torus and shaping the uniform semi-opening angle. 

We thank K. Iwasawa and R. Antomicci for helpful com- 
ments and suggestions that improved the paper. TM is 
a JSPS Research Fellow. This work was financially sup- 
ported in part by Grant-in- Aids for the Scientific Research 
(Nos. 07044054, 10044052, and 10304013) of the Japanese 
Ministry of Education, Science, Sports, and Culture. 
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Fig. 6. — Upper panel: relationships between R{L, 25) and the viewing angle for six dusty torus models given in Table 4. Lower panel: 
distributions of the observed R{L, 25) values. Galaxies shown by open triangles are likely to suffer from contamination and are not used in 
our analysis (see text and Table 1). 
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